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I. INTRODUCTION 

Recently holographic dual models of QCD were pro- 
posed [H, H, H, 13, M H, 01 and applied for investiga- 
tions to the hadronic physics. This approach is based on 
the AdS/CFT correspondence [8j], and is aiming to con- 
struct dual models of QCD in 5-dimensional Anti-de Sit- 
ter (AdS) space. Invented models incorporate important 
features of QCD, as confinement and chiral symmetry 
breaking, and allow one to calculate numerous hadronic 
properties including masses, decay constants, couplings 
of various mesons. 

The holographic models of QCD were employed for in- 
vestigation of mesons electromagnetic form factors (FFs) , 
as well. Calculations of the pion and p-meson FFs in 
Refs. d, [l^, [ll|, [l^, [l^ were carried out in the frame- 
work of two models, namely in hard- wall and soft- wall 
ones depending on the procedure to impose the infrared 
cutoff on the fifth (holographic) dimension of the AdS 
space. 

It is remarkable, that the holographic QCD (HQCD) 
predicts also mesons' valence wave functions thus pro- 
viding analytic approximation to QCD p^ . By this 
way, the HQCD prediction for the pion distribution am- 
plitude (DA) ifhoi{u) was derived. It differs from the 
pion asymptotic DA (pasyiu) found from the perturba- 
tive QCD (PQCD) evolution [TJ], and due to the broader 
shape should increase the magnitude of the leading twist 
QCD prediction for the pion electromagnetic FF. 

In the present work the DA iphoiiu) is used to com- 
pute the pion spacelike electromagnetic form factor in the 
context of the QCD light-cone sum rule (LCSR) method 
with twist-6 accuracy by taking into account 0{as) or- 
der correction to the leading twist term. The twist-4 DA 
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of the pion, necessary to determine the twist-4 contribu- 
tion to the FF, is derived in the context of the renor- 
malon method. The twist-4 contribution is also modeled 
employing the asymptotic form of the pion twist-4 DA 
following from the conformal expansion. 

The paper is structured as follows: In Sec. |TT] general 
expressions for the form factor F-^{Q'^) in LCSR method 
are presented. Section IIIII is devoted to calculation of 
the pion twist-4 DA employing the renormalon approach. 
In Sec. IIVI we present results of numerical computations 
and compare them with the experimental data and pre- 
dictions of the holographic QCD. Section |V] contains our 
conclusions. 



II. THE PION ELECTROMAGNETIC FORM 
FACTOR IN THE QCD LCSR METHOD 

It is known that the QCD LCSR method is one of the 
powerful tools to evaluate nonperturbative components 
of exclusive quantities ■ The LCSR expression for the 
pion electromagnetic FF was derived in Refs. [TgI. [TtI. [lH : 
with the twist-6 accuracy it has the following form 

(Q2) ^ (g2) + (Q2) + ^(4) (q2) ^ ^(6) (q2) ^ 

(2.1) 

where — —q^, q being the four-momentum of the 



virtual photon in the process 7*7r 



*-rr± 



In Eq. 



(HH) Fi"\Q'^) is the twist-n contribution to the FF; 
^^2,Qb)^q2-j jg (3(Q,g) order correction to the twist-2 part. 

The leading twist (twist-2) light-cone sum rule for 
Ftt{Q'^) at the leading order is given by the expression 



where 
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In Eqs. jMl) and jlSl) (p'^'^\u, fi%) is the pion leading 
twist DA; So is the duahty interval, is the Borel vari- 
able and u = 1 — u. 

The 0{as) order correction to the leading twist term 
was obtained in Ref. 17| • Then the whole twist-2 contri- 
bution to the form factor 

F(2)(g2) = Fi2)(g2) + 42,o.)(Q2) 

can be written as the sum of the soft and hard parts 

r-l 



(2.5) 



Here 



(4) / 2\ - (4) ^ 2 N 



(2.6) 

and ¥^2^' (u, is the pion two-particle twist-4 DA. 

The factorizable twist-6 contribution to the LCSR was 
computed in Ref. [l3| in terms of the quark condensate 
density 



Q{u^-u)fI^I^{u,M\s^) 



(2.4) 



The soft part of the LCSR contains entirely the leading 
order twist-2 contribution and some piece of the next- 
to-leading order correction to the twist-2 term. Stated 
differently, the expression presented in Eq. (|2.2|) is the 
pure soft contribution to the FF. The twist-2 term (|2.4p 
is linear in the pion DA. Its hard component at high mo- 
mentum transfers 3> sq leads to the PQCD prediction 
for the FF quadratic in DA For example, choos- 

ing as <p^^^(w, /i|i) the pion asymptotic DA ipasyiu), from 
the hard part of the LCSR one regains the well-known 
asymptotic prediction of the PQCD [3, [Toj 



where Cf = 4/3 is the color factor. 

In the QCD LCSR method for the factorization and 
renormalization scales the following values should be ac- 
cepted: 



(2.8) 



Equations (|2.ip - (|2.7p supplemented by the prescription 
(|2.8[) form a basis for investigation of the pion electro- 
magnetic FF in the QCD LCSR method. 



with = 0.132 GeV being the pion decay constant. 

Further details of calculations and explicit expressions 
for F^{u, M^, So) and F^^J^j^{u, M^, sq) can be found in 
Ref. O. 

The twist-4 term Pj.'^^Q^) is determined as (isj 



du — — exp 

uM2 
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III. THE PION DISTRIBUTION AMPLITUDES 



The light-cone two-particle distribution amplitudes of 
the pion are defined through the light-cone expansion of 
the matrix element 



(0|d(a;2)7.75[x2,a:i]u(xi)|7r+(p)) =z/,p, / due-^f^^'^^P^^ U^\u, fi^) + ^\^^\u, fij,) + 0{A^) 
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(3.1) 



where (pf \u, ^|,), (p'^\u, fij^) are two-particle twist-4 the leading twist pion DA we get [T 
DAs, A = xi — X2, and we use the notation [a;2,a:^i] for 
the Wilson line connecting the points Xi and X2- 



The standard method to handle meson DAs is model- 
ing them employing the conformal expansion. Then for 
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Here (pasy{u) is the PQCD asymptotic DA of the pion 

and Cf/'^{^) are the Gegenbauer polynomials. The func- 
tions bn{ii%) determine the evolution of ip^'^^{u, on 
the factorization scale [3| , 



2^ ^6° 



M) 



In = Cp 



n+1 



(n+ l)(n + 2) 



(3.3) 



In the above 7„ are the anomalous dimensions, /Iq is the 
normalization scale, and fej^ = 6„(/ig). The expansion 



over the conformal spin can also be performed for the 
higher twist DAs (see, for example, Refs. pol [2]|). 

An alternative way to find the higher twist DAs is the 
renormalon approach [2^ , [23j . The renormalon approach 
employs the assumption that the infrared rcnormalons in 
the leading twist coefficient functions should cancel the 
ultraviolet renormalons in the matrix elements of twist-4 
operators in a relevant operator product expansion. Such 
cancellation was proved by explicit calculations in the 
case of the simple exclusive amplitude involving pseu- 
doscalar and vector mesons It turned out that this 
is enough to find the full set of two- and three-particle 
twist-4 DAs of pseudoscalar and vector mesons in terms 
of their leading twist DAs. Higher twist DAs of some of 
the mesons were computed using the renormalon tech- 
nique in the papers 0, [H, . 

In the renormalon-based model the pion twist-4 DAs 
are given by the formulas [1^ 



f't^ {U,fi%)=^-l^ dM^(2) I ^ + _ ^) 1^ _ ^) 
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[v < m)| , 
(3.4) 

(3.5) 



As is seen the renormalon model for the twist-4 DAs 
depends only on one free parameter 5"^. It is related to 
the matrix element of the local operator 



djiyigG^pU 



In Eq. p.7p we use the latest available estimation for 
<5o l25i]. In what follows we also do not show explicitly 
a dependence of the parameter 5"^ on the factorization 



scale 



2^^8C^./3/3o 



(3.6) 



and at /ig = 1 GeV^ was estimated from various 2-point 



QCD sum rules [25|, |2 



(5^(^^) = 5t,= 0.18 ± 0.06 GeV^ 



In the case of the pion leading twist DA with two non- 
asymptotic terms {b2(lJ-'^), b4{pp) ^ 0, and 6„ = 0, n > 
4) these higher twist distributions were calculated in Ref. 
[23 | (see Erratum in Ref. [1^). Let us rewrite the twist-4 
DA ip''2\u, fXp) obtained in Ref. 24|, and solely relevant 



(3.7) for our present studies, in the compact form 
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-1-364 (/^p) ( Sm^ luM + Sm^ lnu + huu -\ m^m^ - 21m^?I^ 
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Then Eq. (j2.6p is given by the expression: 

04 (m, /i^) = 2(5^ {— 4mu — 4mm In M -I- 2m(2u — 1) Inu +462 (/i^) [— IImu + ibu^u — 407i'^M — 6Mulnu + 3m(2u — 1) Ium] 
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-&4(mI) [-137wIt + 917u^M- 3073u^m +4347m%- 2268u^It - 60uulntl + 30u(2u - l)lnu]} . 
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(3.9) 



The distribution amplitude of the pion obtained within 
the HQCD 



(3.10) 



in the renormalon approach leads to the following twist-4 
DA 
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As a result, for the twist-4 function (f>'l"'-{u), we get 
4,52 



3tt 



u(l + 8u)x - + (3 - llii + 8u^)^ 

u \l u 



16uu arctan w — 

u y u 



(3.12) 



The leading twist distributions (pasy{u) and iphoiiu) 
and the corresponding twist-4 DAs ip'^'^ (u) are depicted 
in Fig. [TJ We see that as compared with ipasy{u) the 
holographic DA (phoi{u) is enhanced in the end-point do- 
mains. The difference between the twist-4 DAs shown 
in Fig. HKb) is mild. As a result, the twist-4 functions 
(f)4{u) and 4)\°''{u) that determine the twist-4 contribu- 
tion to the form factor are close to each other (Fig.[2ja)). 
Nevertheless, we should note that the function 4)\°\u) is 
enhanced in the end-point region u 
demonstrate the variation of the form of 4)\°^ 



1. In Fig.lHb) we 

(u) depend- 
ing on the chosen value of the parameter 5q . 

The twist-4 DA of the pion can be modeled in the 
context of the conformal expansion as well. The lowest 
conformal spin, i.e. the asymptotic form of the twist-4 
DA (^2''^(u) was employed in Ref. [l^ for calculation of 
the twist-4 function (|2.6p with the following result 



(3.13) 



90 

Cj3l'y(u) = yfe[l-u(7-8u)]. 



We use the function 04*^ (u) for calculation of the asymp- 
totic twist-4 term. 



IV. NUMERICAL RESULTS 

In order to perform numerical computations we should 
fix various parameters appearing in the relevant expres- 
sions. Namely, we take the Borel parameter equal to 



A/^ — 1 GeV^ and accept for the factorization and renor- 
malization scales the values presented in Eq. (|2.8p . For 
the QCD coupling as(/ip) the two-loop expression with 
A3 = 0.34 GeV is used. The value of the duahty pa- 
rameter So = 0.7 GeV'^ is borrowed from the QCD sum 
rule for the correlator of two u'jfj.jsd currents [29j. The 
normalization scale is set equal to /Xg = 1 GeV^. 

Results of numerical computations are depicted in 
Figs. Hi As is seen (Fig. (SJa)) the leading twist LCSR 
contribution found employing iphoi{u), due to the broader 
shape of the HQCD distribution amplitude, is larger than 
the prediction obtained by means of ipasy{u). The twist- 
4 terms computed using the twist-4 function ^4°' (u) and 
the asymptotic version of Eq. p.9p do not differ consid- 
erably from each other: one observes some deviation of 
the solid curve from the result of the standard QCD, i.e. 
from the dashed line. Here, the features of the different 
twist-4 terms that have interesting consequences should 
be emphasized. Thus, the asymptotic twist-4 term (the 
dot-dashed line) in the region of low momentum trans- 
fers demonstrates more rapid increasing on than other 
twist-4 terms: in this region among the twist-4 contribu- 
tions the HQCD prediction is a flattest and lowest one. 
Contrary, for the momentum transfers > 8 GeV^ 
the HQCD twist-4 term overtakes other twist-4 contribu- 
tions. What is important, in the whole region of consid- 
ering momentum transfers only the HQCD twist-4 term 
is monotonically increasing function of Q^. All twist-4 
terms are sensitive to the choice of the parameter Sq. In 
Fig. [3Kb), as an example, we plot the twist-2 and -4 terms 
generated by the HQCD distributions and by shaded area 
show sensitivity of the twist-4 term to 5^ . 

The pion scaled electromagnetic FF Q'^F.^{Q^) com- 
puted employing the HQCD distribution amplitude and 
twist-4 function 04°' (w) is depicted in Fig. UK a). It is 
clear that the enhancement in the leading twist term is 
not sufficient to describe the existing data on F^(Q^); the 
corresponding LCSR curve runs below the data points. 
If we replace the holographic twist-4 term by its asymp- 
totic counterpart keeping the remaining ones unchanged, 
then the final expression describes the experimental re- 
sults. Nevertheless such agreement can be achieved only 
for the large values of the parameter 5^ > 0.2. From this 
analysis it is legitimate to conclude that the enhancement 
generated by the shape of the function iphoi (u) in the rel- 
evant terms of Eq. (|2.1|) and the asymptotic twist-4 term 
is enough to explain the data on F^(Q^) in the LCSR 
treatment. 

The pion electromagnetic form factor was analyzed 
within the LCSR method in Ref. [24|, where in calcu- 
lations model DAs with two nonasymptotic terms and 
renormalon inspired twist-4 DAs were used, and from 
comparison with the data, constraints on the input pa- 
rameters 62 and 64 were extracted. In the fitting proce- 
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FIG. 2: (a) The twist-4 function <j)4,{u) obtained in the 
renormalon-based model using the PQCD asymptotic DA 
(the dashed line) and the distribution amplitude of the HQCD 
(the solid line); (b) the dependence of the twist-4 function 
04°' (w) on the input parameter ^g- 



dure in Ref. [1^ both the old experimental results [s^l 
and new ones from the Jefferson Lab. Collaboration 
[sij l were used. It is worth noting that the pion form 
factor is not measured directly and real measurements of 
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FIG. 3: (a) The leading and twist-4 contributions to the 
scaled pion form factor. In the calculations the PQCD DAs 
(the dashed lines) and DAs of the HQCD (the solid lines) 
are used. For comparison the asymptotic twist-4 term (the 
dot-dashed line) is also shown. All twist-4 terms are obtained 
using 5q = 0.2. (b) The leading and twist-4 contributions 
to the FF found employing the functions iphoi{u) and </)J°'(it). 
The shaded area demonstrates the allowed values of the twist- 
4 term obtained by varying the parameter 5q within the limits 
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FIG. 4: (a) The pion scaled FF as a function of . It is com- 
puted using the HQCD distribution amplitude iphoi{u) and 
the renormalon-based twist-4 function <j)\°''{u). The shaded 
area shows allowed limits of the form factor. For the cen- 
tral solid line Sl = 0.2. (b) The pion FF obtained using the 
asymptotic twist-4 term and the holographic ones for the re- 
maining contributions from Eq. (|2.ip . For the central solid 
line 5q = 0.2. Concerning the data - see discussion in the text 



the process 7*p 



n require an extrapolation of the 



off-shell pion to the mass-shell. Because this extrapo- 
lation becomes increasingly problematic as increases 
and at high momentum transfers the old results suffer 
from large uncertainties, in the present work we remove 
them from consideration. We also correct the data from 
the i^TT Collaboration reported in Ref. [13], and add two 
new points extracted by Ft^2 Collaboration [s^: in Figs, 
m [5] and [5] the rectangle and triangle symbols denote the 
F-jj and ¥^^2 data, respectively. Strictly speaking, the 
LCSR method is applicable in the domain of momen- 
tum transfers > 1 GeV^. Nevertheless, we extend 
our numerical calculations to the region < 1 GeV^ to 
reveal properties of the leading and twist-4 terms at low 
momentum transfers. The data below this border (the 
solid circles) are included into the figures for illustrative 
purposes and are irrelevant for our present discussions. 

We have repeated calculations of Ref. [23| using the 
pion model DAs with two nonasymptotic terms. One 
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FIG. 5: The pion scaled FF as a function of . The cen- 
tral solid line is the LCSR prediction obtained employing the 
PQCD leading twist DA with bl = 0.25, h\ = Q and corre- 
sponding renormalon-based twist-4 function cj>4{u,Q^). The 
shaded area is obtained by varying Sq within the allowed lim- 
its at fixed 6° = 0.25, 6° = 0; for the central hne Sq = 0.2. 
For comparison the LCSR result found using ipaay(u) is also 
shown (the dashed line, 5q — 0.2). 



of these FFs, that corresponds to values of the input 
parameters 5° = 0.25, 64 = and (5§ = 0.2, is plotted 
in Fig. [5l The achieved agreement with new data of the 
JLAB experiment is encouraging. It is easy to see that 
excluding of the unprecise data from analysis and adding 
new ones results in a shift of the parameter towards 
larger values than in Ref. [24] . The similar curves can be 
obtained in the more general case 63 7^ 0, 64 7^ as well. 

In the framework of the holographic QCD the 
electromagnetic FF has been calculated in Refs. 
[T3| . The soft- wall model expression for the form factor 
is especially simple, and is given by the monopole form 



4fc2 



4fc2 + Q2 



(4.1) 



where k = 0.375 GeV. 

In Fig. [S] we compare the LCSR predictions found using 
the functions iphoi{u) and (f)\°^[u) (the lower solid hne) 
and ifhoi{u) and (j)2''^{u) (the dashed line) with the pre- 
diction of the soft-wall HQCD (the upper solid line) . Let 
us note that the dashed line has been obtained choosing 
the maximal allowed value of — 0.24. 

The first LCSR result with the HQCD distribution am- 
plitude and renormalon inspired twist-4 function in the 
range of the momentum transfers 2 < < 10 GeV^ 
demonstrates the scaling in accordance with Q^Ft^{Q'^) ~ 
const law with a weak sign of the logarithmic running 
of the QCD coupling. This feature can be explained 
as a consequence of the dependence of the correspond- 
ing twist-4 term on due to the enhancement of the 
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FIG. 6: The pion scaled FF as a function of Q^. The lower 
solid line is prediction of the LCSR method with the distribu- 
tion amplitude of the HQCD. The dashed line is the same the 
HQCD result, but the HQCD twist-4 term replaced by the 
asymptotic twist-4 one. The upper solid line is the soft-wall 
holographic QCD prediction from Ref. p^ . 

(y94°'(u) at end-point regions. The second LCSR predic- 
tion computed also in this work differs from the first one 
only in the choice of the twist-4 term. It is evident that 
the dashed line shows more rapid fall-off with than 
the lower solid line, though due to the choice Sq = 0.24 is 
in agreement with the data. Such rapid falling is typical 
for other LCSR curves obtained using the PQCD DAs 
(Fig.©. 

In the LCSR method the soft component of the pion 
FF is modeled by the leading order twist-2 term and 
some piece of the 0{as) order correction to the twist- 
2 term. The leading order higher twist terms are soft 
contributions to the pion FF as well. There are various 
approaches to model the soft component of the FF. For 
example, in Ref. [sj the pion FF was evaluated apply- 
ing the standard hard-scattering approach and the QCD 
running coupling method. It was demonstrated that the 
Borel resumnied expression of the form factor contains 
both the soft and hard components, and in the asymp- 
totic limit Q'^ ^ oo the standard PQCD result can be 
recovered. Alternatively, in Ref. [s^ soft and hard contri- 
butions to the pion electromagnetic FF were calculated 
within the light-front quark model and comparison with 
the AdS/CFT prediction was performed. The authors 
reported on excellent agreement of their result with the 
HQCD prediction. 

Comparison of our results with the soft-wall holo- 
graphic QCD prediction reveals the discrepancy between 
them. First of all, the prediction of the holographic QCD 
overestimates the JLAB data: a situation with a hard- 
wall result is even worse fT2|. Second, the HQCD predic- 
tion obeys the scaling Q^Ft^IQ"^) const with n ^ 2. 
From our point of view, emphasized also in Ref. (ssj . for a 
reliable analysis one should evolve the pion DA to a scale 



of the considering process, which in our case is equal to 
a few GeV^. This procedure inevitably demands to in- 
clude nonasymptotic terms into the holographic DA of 
the pion. In other words, for credible phenomenological 
applications the holographic QCD version of Eq. (|3.2p . 
where the function iph,oi{u) is only the first term in a rel- 
evant expansion, is required. Higher order corrections 
neglected in present investigations of the AdS/CFT cor- 
respondence may improve a situation, too [s^ . 

For successful applications of the pion HQCD distribu- 
tion amplitude within the framework of traditional meth- 
ods of the perturbative QCD, one should prove that this 
DA describes also other exclusive and semi-inclusive pro- 
cesses involving the pion. The simplest such process is 
the pion electromagnetic transition form factor Ft^^{Q'^). 
In the context of the LCSR method it was computed in 
Refs. [sgI . [37} : the result was derived with twist-4 ac- 
curacy including the next-to-leading order correction to 
the twist-2 term. In the context of the LCSR and renor- 
malon methods this form factor was analyzed in Refs. 
(28I [sst . Preliminary calculations demonstrate that the 
twist-2 term in the LCSR expression for F.^j{Q'^) found 
using fhoiiu) exceeds the measured experimental data 
[3^. But contributions of the next-to- leading order and 
twist-4 terms to the LCSR are negative and reduce the 
magnitude of the leading twist contribution (see, for ex- 
ample, [28j): hence an agreement with the data may 
be achieved. A complete analysis of F-^^iCp) requires 
to compute using iPhni(u) the renormalon-based twist-4 
function (u, Q^) [H, [si] that determines the twist-4 
contribution to the transition FF. This problem requires 
separate investigations. 



V. CONCLUSIONS 

In the present work we have calculated the pion space- 
like electromagnetic form factor Ft^{Q^) in the QCD 
LCSR framework. We have used the distribution ampli- 
tude of the pion derived in the holographic QCD, which 
has a broader shape than the PQCD asymptotic DA. The 
twist-4 DA has been obtained applying the renormalon 
inspired model of Ref. |23|. The prediction of the LCSR 
method obtained with the twist-6 accuracy by means of 
the functions (phoi{u) and (j)'l°^{u) lies below the data, 
and does not describe the experimental situation. The 
agreement with the data can be achieved provided the 
HQCD twist-4 term is replaced by the asymptotic one 
and (5q > 0.2 is chosen. 

We have also computed F^(Q^) utihzing the pion 
PQCD leading twist DAs with two nonasymptotic terms 
and corresponding renormalon inspired twist-4 DAs. It 
has been proved that in this case in the region > 
1 GeV^ the QCD LCSR method explains the new exper- 
imental situation emerged due to the data of the JLAB 
experiment. 

Comparison with results of the holographic QCD [ll|, 
[T2| . [Tst has revealed an interesting tendency: the holo- 
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graphic QCD predictions for Q^F-^{Q'^) obtained within 
the both hard- and soft-wall models are increasing func- 
tions of Q^, whereas in the LCSR treatment at large 
momentum transfers one finds an opposite picture. It 
is worth noting that the scaling of the form factor cal- 
culated employing the functions fhoiiu) and tf\°\u) is 
very promising. Unfortunately, the corresponding curve 
runs below the data. In order to enhance the magni- 



tude of Ft^{Q'^) and reach an agreement with the data 
additional contributions are needed. It seems to us that 
such contributions may appear due to the holographic 
QCD counterparts of the nonasymptotic terms in the 
pion PQCD leading twist DA (|3.2[) . If extracted, these 
terms may improve scaling properties and normalization 
of the holographic QCD predictions themselves. 
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